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ABSTRACT 

A  geometric  model  has  been  developed  for  analyzing  observed  Beaufort/Chukchi  Sea  temperature  and  salinity 
profiles  based  on  their  geometric  features.  It  contains  four  major  components:  (a)  a  nonhomogeneous  mesh 
difference  scheme  for  second  derivatives  d2T{zk)/dz 2  and  d2S(zk)ldz2  at  each  data  point  (z*,  Tk,  Sk );  (b)  a  set  of 
key  depths  (depths  of  mixed  layer,  upper  and  lower  thermoclines  and  haloclines,  and  deep  layers)  determined 
mostly  by  the  geometric  features  of  profiles,  such  as  depths  of  minimum  and  maximum  of  d2T(zk)/dz2  and 
d2S(zk)/dz2  for  shelf  water  and  upper  portion  deep  water;  (c)  modeled  profiles  T(zk)  and  S(zk),  which  are  linearly 
interpolated  from  observed  T,  S  values  at  the  key  depths;  and  (d)  statistical  tests  for  accepting  [T^z*),  5(zA.)]  as 
a  good  representative  of  ( zk ,  Tk,  Sk ).  The  output  of  this  geometric  model  is  a  set  of  major  characteristics  of  each 
profile:  sea  surface  temperature  (salinity),  mixed  layer  depth,  upper-  and  lower-thermocline  (halocline)  depths, 
thermocline  (halocline)  temperature  (salinity)  gradient,  and  deep  layer  stratification.  Analyzing  3562  Beaufort/ 
Chukchi  Sea  historical  (1970-93)  temperature  and  salinity  profiles  from  the  Naval  Oceanographic  Office’s  Master 
Oceanographic  Observation  Data  Set  by  this  geometric  model,  the  Beaufort/Chukchi  Sea  thermohaline  field 
reveals  a  dual  structure:  one  layer  (vertically  uniform)  on  the  continental  shelf  (depth  <  150  m)  during  winter 
and  multiple  layers  (mixed  layer,  thermocline/halocline,  sublayer)  during  summer.  Strong  seasonal  variations 
were  also  found  in  mixed  layer  depth,  upper  thermocline/halocline  depth,  and  strength. 


1.  Introduction 

Determination  of  thermohaline  structure  (mixed  layer, 
entrainment  zone,  thermocline,  halocline)  from  ob¬ 
served  temperature  and  salinity  profiles  is  quite  impor¬ 
tant  in  the  Arctic  Ocean  for  several  reasons.  First,  the 
heat  balance  of  the  Arctic  Ocean  depends  on  the  features 
of  the  mixed  layer,  entrainment  zone,  and  thermocline. 
The  mixed  layer  deepens  by  entrainment  of  water  from 
the  ocean  below  (Chu  and  Garwood  1991).  The  heat  in 
the  deep  water  is  then  transported  to  the  surface.  Second, 
the  mixed  layer  acts  as  a  buffer  by  storing  solar  heat 
input  during  the  summer  and  releasing  it  back  to  the  ice 
throughout  the  fall  and  early  winter  (Maykut  and 
McPhee  1995). 

Teague  et  al.  (1990)  developed  a  curve-fitting  method 
to  build  upper-ocean  (0-400  m)  temperature  and  salinity 
climatology  for  the  Generalized  Digital  Environmental 
Model  (GDEM).  The  parameters  they  used  are  purely 
mathematical.  To  get  physically  meaningful  parameters. 
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Chu  (1995)  and  Chu  et  al.  (1997)  proposed  a  parametric 
model  based  on  a  single  type  profile  (mixed  layer,  en¬ 
trainment  zone,  thermocline,  deep  layer)  and  an  iteration 
scheme.  The  model  consists  of  three  major  components: 
(a)  a  first-guess  model  profile,  (b)  high-resolution  pro¬ 
files  interpolated  from  observed  profiles,  and  (c)  fitting 
of  high-resolution  profiles  to  the  parametric  model  pro¬ 
file  through  iteration.  The  output  is  a  set  of  major  char¬ 
acteristics  of  each  profile:  sea  surface  temperature, 
mixed  layer  depth,  thermocline  depth,  thermocline  tem¬ 
perature  gradient,  and  deep  layer  stratification. 

The  Chukchi  Sea  is  a  shallow  sea  with  a  mean  depth 
of  40-50  m,  having  gentle  knolls  and  several  troughs 
that  are  shallow  but  with  a  relief  that  is  a  substantial 
fraction  of  the  mean  depth.  The  Beaufort  Sea  is  a  semi- 
enclosed  basin  with  a  quite  narrow  continental  shelf 
(30-80  km)  and  a  deep  submarine  canyon,  the  Mac¬ 
kenzie  Canyon.  The  deepest  water  is  confined  to  the 
bowl-shaped  Mackenzie  Canyon,  which  is  enclosed  by 
the  3000-m  isobath.  The  gradients  in  slope  across  the 
bottom  are  very  large.  The  150-m  isobath  approximately 
characterizes  the  Beaufort  Sea  shelf  break  (Fig.  1).  The 
shallow  continental  shelf  waters  are  readily  affected  by 
seasonally  varying  atmospheric  conditions,  such  as 
heating,  cooling,  wind  stress,  and  ice  formation  and 
melting.  The  seasonal  variation  of  the  water  masses  is 
remarkably  large.  Thus,  we  have  two  different  types  of 
profiles:  (a)  shelf  profiles  (water  depth  <  150  m)  and 
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Fig.  1.  The  Beaufort/Chukchi  Sea  bathymetry.  The  data  were  obtained  from  the  Naval  Oceanographic  Office  DBDB5  world  bathymetry 

database.  Numbers  show  the  depth  in  100  m. 


(b)  deep  water  profiles  (water  depth  >150  m).  Among 
the  historical  (1970-93)  temperature  and  salinity  pro¬ 
files  (3562)  obtained  from  the  Naval  Oceanographic  Of¬ 
fice,  3384  profiles  are  shelf  profiles  and  178  profiles  are 
deep  water  profiles. 

The  usefulness  of  the  parameter  model  (Chu  1995; 
Chu  et  al.  1997)  is  based  on  the  single-type  approxi¬ 
mation  on  profiles.  If  T,  S  profiles  follow  the  assumed 
pattern,  the  parametric  model  has  a  good  capability  to 
grasp  major  features  of  observed  profiles.  If  T,  S  profiles 
demonstrate  multiple  types,  it  is  difficult  to  use  the  para¬ 
metric  model  since  the  first-guess  model  profile  is  hard 
to  determine.  Thus,  in  this  study  we  developed  a  new 
model  to  process  multiple-type  T,  S  profiles  in  the  Chuk¬ 
chi/Beaufort  Sea  based  on  the  geometric  features  of  the 
T,  S  profiles,  such  as  the  depths  of  the  minimum  and 
maximum  second  derivatives  d2T/dz2,  d2S/dz2.  With  us¬ 
ing  the  observed  T,  S  values  at  these  depths  (called  key 
depths)  and  omitting  any  first-guess  modeled  profile, 
the  geometric  model  transforms  each  profile  into  a  set 
of  characteristic  parameters  such  as  sea  surface  tem¬ 
perature  (SST),  sea  surface  salinity  (SSS),  mixed  layer 
depth  (MLD),  temperature  and  salinity  jumps  across  the 
entrainment  zone  (AT,  AS),  upper-  and  lower-thermo- 
cline  (halocline)  depths  and  strengths,  and  the  deep  layer 
stratification. 

2.  Master  Oceanographic  Observation  Data  Set 
(MOODS) 

The  MOODS  is  a  set  of  observed  ocean  data  world¬ 
wide  consisting  of  (a)  temperature-only  profiles,  (b) 


both  temperature  and  salinity  profiles,  (c)  sound-speed 
profiles,  and  (d)  surface  temperatures  from  drifting 
buoys.  It  includes  observations  from  various  instru¬ 
ments  such  as  conductivity-temperature-depth  (CTD), 
salinity-temperature-depth  (STD),  low-resolution  STD, 
expendable  bathythermograph,  and  so  on.  These  mea¬ 
surements  are,  in  general,  irregular  in  time  and  space. 
The  data  are  compiled,  archived,  and  updated  by  the 
Naval  Oceanographic  Office.  In  this  study,  we  analyze 
temperature  and  salinity  profiles  acquired  from  CTD/ 
STD  and  bottle  data.  Our  study  domain  includes  the 
area  65°-75°N,  120°-180°W;  the  dataset  within  this  re¬ 
gion  during  1970-93  consisted  of  3562  profiles  after 
rejecting  certain  data  during  quality  control.  These  pri¬ 
mary  editing  procedures  included  removal  of  profiles 
with  obviously  erroneous  location,  profiles  with  large 
spikes,  poor  vertical  resolution,  and  profiles  displaying 
features  that  do  not  match  the  characteristics  of  sur¬ 
rounding  profiles. 

The  main  limitation  of  the  MOODS  data  is  its  irreg¬ 
ular  distribution  in  time  and  space.  Certain  periods  and 
areas  are  well  sampled,  while  others  lack  enough  ob¬ 
servations  to  gain  any  meaningful  insights.  Vertical  res¬ 
olution  and  data  quality  are  also  highly  variable,  de¬ 
pending  to  a  large  extent  on  instrument  type  and  sam¬ 
pling  expertise.  There  is  a  data  sparse  area  north  of  73°N 
(Fig.  2a).  The  periods  of  1971-72,  1977,  1979,  and 
1985-89  are  found  to  have  a  relatively  large  number 
profiles  averaging  around  200  profiles  per  year  (Fig.  2b). 
Most  profiles  were  observed  during  the  summer  season 
(Fig.  2c).  August  logs  the  most  observations  (1253), 
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Fig.  2.  MOODS  data  for  this  study:  (a)  distribution  of  stations,  (b)  interannual  variation  (1970-93)  of  number  of  profiles,  and  (c)  season¬ 
al  variation  of  number  of  profiles. 


while  no  observations  have  been  recorded  in  the  De¬ 
cember  and  January  dataset. 

3.  Seasonal  variation  of  temperature  and  salinity 

profiles 

a.  Basic  features 

The  surface  radiative  flux  strongly  affects  the  Beau¬ 
fort/Chukchi  Sea  shelf  thermohaline  structure.  During 
the  winter  (November-April),  long  periods  of  darkness 
together  with  low  solar  elevation  give  rise  to  a  prolonged 
period  of  radiative  loss  from  the  surface.  The  radiative 
cooling  at  the  surface  destabilizes  the  upper  layer 
through  strong  upward  heat  flux  and  salt  rejection  by 
ice  freezing,  and  causes  the  formation  of  a  deep  mixed 
layer,  which  on  the  shallow  shelves  might  reach  the 
bottom  resulting  in  an  isothermal/isohaline  structure 
(Figs.  3a, b).  Different  lengths  of  these  profiles  in  the 
vertical  are  caused  by  the  different  water  depths  where 
the  observations  were  taken  (Fig.  1).  During  the  summer 
(August),  long  daylight  hours  together  with  relatively 
high  solar  elevation  give  rise  to  a  period  of  radiative 


deposition  to  the  surface.  The  radiative  warming  at  the 
surface  stabilizes  the  upper  layer  through  downward 
heat  flux  and  freshwater  influx  by  ice  melting,  causing 
the  mixed  layer  to  shoal  and  forming  a  multilayered 
structure  (i.e.,  a  mixed  layer,  upper  and  lower  thermo- 
clines  and  haloclines,  and  a  deep  layer),  as  shown  in 
Figs.  3d,e.  The  number  and  spatial  distribution  of  ob¬ 
servations  are  much  greater  in  summer  than  in  winter 
(Figs.  3c,f). 

In  the  deep  water  off  the  shelf  the  near-surface  waters 
of  the  Beaufort  Sea  also  experience  seasonal  variations. 
During  the  winter  (November-April),  surface  cooling 
causes  the  formation  of  a  deep  thermal  mixed  layer  (Fig. 
4a).  However,  the  surface  salt  flux  caused  by  local  ice 
freezing  generates  a  relatively  shallow  salinity  mixed 
layer  (Fig.  4b).  Below  both  the  thermal  and  salinity 
mixed  layers  there  exists  a  lower  thermocline  and  halo- 
cline,  appearing  at  160-300-m  depth  (Figs.  4a, b).  Dur¬ 
ing  the  summer  (August),  surface  warming  and  asso¬ 
ciated  ice  melting  increase  the  SST  (a  maximum  value 
near  8°C),  decrease  the  SSS  (a  minimum  value  near  20 
psu),  and  cause  both  the  thermal  and  salinity  mixed 
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Fig.  3.  Seasonal  variation  of  T,  S  profiles  over  the  Beaufort/Chukchi  Sea  shelf  water:  (a)  winter  (Nov- Apr)  T  profiles,  (b)  winter  S  profiles, 
(c)  distribution  of  winter  stations,  (d)  summer  (Aug)  T  profiles,  (e)  summer  S  profiles,  and  (f)  distribution  of  summer  stations.  The  shaded 
areas  indicate  shelf  region  (water  depth  <150  m). 


layers  to  shoal  (Figs.  4d,e).  We  also  notice  that  both 
winter  and  summer  stations  (Figs.  4c, f)  do  not  extend 
far  from  the  shelf  break,  and  the  intermediate  waters 
(below  150-m  depth)  do  not  exhibit  a  seasonal  variation. 

b.  Classification  of  T,  S  profiles 

The  historical  Beaufort/Chukchi  Sea  temperature  and 
salinity  profiles  (Figs.  3  and  4)  demonstrate  the  exis¬ 
tence  of  several  basic  profile  shapes,  which  we  term 
deep  mixing  (DM),  shallow  mixing  (SM),  and  advection 
(AD). 

1)  Deep  mixing  (DM)-type  T,  S  profiles 

The  DM-type  T,  S  profiles  are  characterized  by  a 
single  well-mixed  layer  in  the  shelf  region  (Fig.  5a)  and 
a  deep  (>150  m)  surface  thermal  mixed  layer  (Fig.  6a) 
and  a  thermocline  in  the  deep  water  region.  This  type 
of  profile  is  caused  by  surface  destabilization,  such  as 
strong  wind  forcing,  surface  cooling,  and  brine  rejection 
due  to  ice  freezing.  They  are  most  prevalent  during 
winter. 


2)  Shallow  mixing  (SM)-type  T,  S  profiles 

The  SM-type  T,  S  profiles  are  characterized  by  a  mul¬ 
tilayered  structure  (a  shallow  mixed  layer,  entrainment 
zone,  thermoclines  and  haloclines,  and  a  sublayer)  in 
both  the  shelf  region  (Figs.  5b, d)  and  the  deep  water 
region  (Figs.  6b, d).  This  profile  type  is  generally  present 
during  the  summer  and  is  caused  by  surface  stabiliza¬ 
tion,  such  as  weak  wind  forcing,  surface  warming,  and 
freshwater  influx  due  to  ice  melting  and  river  runoff, 
and  therefore  usually  has  a  warm  SST  and  a  low  SSS. 
Beneath  the  mixed  layer  the  thermocline/halocline  pro¬ 
file  is  complex,  reflecting  the  seasonal  adjustment.  Thus, 
SM  profiles  feature  a  monotonic  decrease  of  temperature 
with  depth  in  the  upper  thermocline  and  a  monotonic 
increase  of  salinity  with  depth  in  the  upper  halocline 
(Figs.  5b,d  and  6b, d). 

3)  Advection  (AD)-type  T  profiles 

The  AD-type  T  profiles  also  have  a  multilayer  struc¬ 
ture  (a  shallow  mixed  layer,  entrainment  zone,  upper 
and  lower  thermoclines,  and  a  sublayer)  in  both  the  shelf 
region  (Fig.  5c)  and  the  deep  water  region  (Fig.  6c). 
This  type  of  profile,  identified  by  a  nose-shaped  curve 
in  the  upper  thermocline  (monotonic  increasing  T  with 
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Fig.  4.  Seasonal  variation  of  T,  S  profiles  over  the  Beaufort/Chukchi  Sea  deep  water:  (a)  winter  (Nov-Apr)  T  profiles,  (b)  winter  S  profiles, 
(c)  distribution  of  winter  stations,  (d)  summer  (Aug)  T  profiles,  (e)  summer  S  profiles,  and  (f)  distribution  of  summer  stations.  The  shaded 
areas  indicate  shelf  region  (water  depth  £  150  m). 


depth  to  a  maximum  value  7m.lx  and  then  monotonic 
decreasing  with  depth),  is  generated  by  a  prominent 
hydrographic  feature  on  the  Beaufort  Sea  shelf:  a  sub¬ 
surface  temperature  maximum,  generally  found  at  about 
20-30-m  depth  in  the  deep  water  and  at  about  10-m 
depth  in  the  shelf  water.  This  temperature  maximum  is 
associated  with  the  eastward-flowing  Bering  Sea  water 
(Coachman  and  Barnes  1961;  Aagaard  1984).  The  warm 
water  that  enters  the  Beaufort  Sea  comes  through  the 
eastern  Bering  Strait  and  follows  the  Alaskan  coast 
around  Point  Barrow.  Mountain  et  al.  (1976)  identified 
this  warm  intrusion  as  the  combination  of  the  two  water 
masses:  warm  (5°-10°C)  and  fresh  (salinity  below  31.5 
psu)  Alaskan  Coastal  Water  (ACW),  and  the  more  saline 
Bering  Sea  Water  (BSW).  The  two  water  masses  mix 
rapidly  with  the  ambient  surface  water  as  they  move 
eastward.  Aagaard  (1984)  estimated  that  ACW  is  not 
clearly  identifiable  east  of  147°-148°W,  and  BSW  east 
of  143°W.  Thus,  in  this  study  we  choose  the  longitude 
of  141°W  as  a  separation  for  the  western  and  eastern 
parts  of  the  Beaufort  Sea  shelf  water. 

c.  Occurrence  of  various  types  on  continental  shelf 

The  Beaufort/Chukchi  Sea  shelf  water  demonstrates 
a  strong  seasonal  variation.  During  winter  (November- 
April),  most  temperature  and  salinity  profiles  are  of  the 


DM  type  (Figs.  3a,b,  and  5a).  During  summer  (August), 
temperature  and  salinity  profiles  (Figs.  3d,e)  indicate 
multilayer  structure,  that  is,  SM-type  T,  S  profiles  and 
AD-type  T  profiles.  To  further  illustrate  the  spatial  var¬ 
iation  of  the  thermohaline  structure,  we  divided  the  shelf 
region  into  three  parts  (Fig.  7):  (a)  southern  Chukchi 
Sea  (65°-70°N,  160°W-180°);  (b)  northern  Chukchi/ 
western  Beaufort  Sea  shelf  (70°-75°N,  141oW-180°); 
and  (c)  eastern  Beaufort  Sea  shelf  (68°-75°N,  120°- 
141°W).  The  relative  number  of  summer  profiles  in  each 
of  the  categories  is  listed  in  Tables  1  and  2. 

In  August,  as  expected,  the  DM  type  of  T  profiles 
exhibit  less  occurrence  (23%  in  total  region,  only  1% 
in  the  eastern  Beaufort  shelf)  than  the  SM  type  of  T 
profiles  (Table  1).  The  number  of  SM-type  profiles  is 
much  more  than  the  AD-type  profiles.  Most  AD-type  T 
profiles  (174)  appear  in  the  northern  Chukchi/western 
Beaufort  shelf,  where  the  Bering  Sea  water  intrusion 
appears  to  flow  over  the  broad  shelf  poleward  into  deep¬ 
er  water. 

In  August,  the  DM  type  of  S  profiles  (Table  2)  is  less 
prevalent  than  the  SM  type  in  the  northern  Chukchi/ 
western  Beaufort  Sea  shelf  (DM  type:  19%;  SM  type: 
81%)  and  in  the  eastern  Beaufort  Sea  shelf  (DM  type: 
8%;  SM  type:  92%).  However,  the  DM  type  dominates 
the  southern  Chukchi  Sea  (DM  type:  77%;  SM  type: 
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(a)  (b) 


T  S 


Fig.  5.  Typical  T,  S  profiles  for  the  shelf  water:  (a)  DM-type  T,  S\  (b)  SM-type  T\  (c)  AD-type  T\  and  (d)  SM-type  S. 


23%).  The  DM-type  occurrence  decreases  from  the 
southern  Chukchi  Sea  to  the  northern  Chukchi/West 
Beaufort  Sea  shelf,  and  further  reduces  toward  the  east¬ 
ern  Beaufort  Sea  shelf. 

d.  Occurrence  of  various  types  in  deep  water 

The  Beaufort  Sea  deep  water  also  experiences  a  sea¬ 
sonal  variation,  but  not  as  strong  as  the  shelf  water.  Most 


winter  (December- April)  historical  (1970-93)  data 
show  DM-type  temperature  profiles  (Table  3).  On  the 
other  hand,  all  three  types  of  T  profiles  are  present  in 
the  August  historical  (1970-93)  temperature  data  (Table 
3).  Among  them,  the  DM  type  accounts  for  45%,  the 
SM  type  for  39%,  and  the  AD  type  for  16%.  Regardless 
of  season,  the  salinity  generally  exhibits  an  SM-type 
profile,  that  is,  nearly  isohaline  with  a  weak,  near-sur- 
face  gradient  sometimes  present. 
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T  and  S  Typical  Profiles  in  Deep  Water 


Fig.  6.  Typical  T,  S  profiles  for  the  deep  water:  (a)  DM-type  T,  S;  (b)  SM-type  T\  (c)  AD-type  T\  and  (d)  SM-type  S. 


4.  Geometric  features  of  T,  S  profiles 

Second  derivatives  d2F(z,)/dz2  (F  is  T  or  S)  are  used 
below  to  describe  the  geometric  features  of  profiles  and 
to  identify  various  types  of  T,  S  profiles  in  the  Beaufort 
and  Chukchi  Seas. 

a.  Nonliomogeneous  mesh  difference  scheme 

For  each  profile  (temperature  or  salinity),  F(zf,  i  = 
1,  2,  .  .  .  ,  N,  we  use  the  following  nonhomogeneous 


mesh  difference  scheme  for  the  second-order  deriva¬ 
tives: 


d2F 

~  1 

( Fi+i  -  F,  F,  -  Ft\ 

dz 2 

Zi  Zi+ 1  —  Zj- 1 

\  1  %i  Zi-  ij 

Here,  i  =  1  refers  to  the  surface,  with  increasing  values 
indicating  downward  extension  of  the  measurement. 
Equation  (1)  shows  that  we  need  two  neighboring  val¬ 
ues,  Fi_l  and  Fi+l,  to  compute  the  second-order  deriv- 
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Longitude(E),  Deplh=0(m) 

Fig.  7.  Division  of  Beaufort/Chukchi  Sea. 


ative  at  z.r  For  i  =  1,  N,  we  use  the  next  point  value, 
that  is, 

d2F 
dz2 

After  the  second-order  difference  is  computed,  we  can 
use  the  following  model  to  determine  the  top  and  bottom 
of  the  upper  thermocline  and  halocline. 

The  mixed  layer  depth  has  a  large  variability,  rang¬ 
ing  from  few  meters  to  tens  of  meters.  To  resolve  the 
mixed  layer  depth  from  T,  S  profiles,  we  have  to  use 
T,  S  values  one  depth  increment  above  and  below  for 
computing  the  second  derivatives.  This  is  due  to  that 
the  finest  vertical  resolution  is  1  m  in  our  dataset  (CTD 
data). 

b.  Features  of  the  shelf  water  profiles 

As  stated  previously,  there  are  three  types  of  T profiles 
(DM,  SM,  and  AD  types)  and  two  types  of  S  profiles 
(DM  and  SM  types)  found  in  the  shelf  water.  The  DM- 
type  T,  S  profiles  reveal  a  single-layer  structure  (Fig. 
5a);  that  is,  the  surface  mixed  layer  extends  to  the  bot¬ 
tom  of  the  shelf  (vertically  uniform). 

The  SM-type  T  profiles  are  characterized  by  a  min¬ 
imum  (maximum)  of  d2T/dz2  at  top  (bottom)  of  the 
thermocline.  Thus,  the  depths  of  the  first  minimum  and 
the  first  maximum  of  d2T/dz2  are  taken  as  the  top,  df, 
and  bottom,  df,  of  the  thermocline  (Fig.  5b),  respec¬ 
tively. 

The  AD-type  T  profiles  are  depicted  by  maxima  of 
d2T/dz2  at  both  top  of  the  upper  thermocline  and  bottom 
of  the  lower  thermocline  and  a  depth,  d^x,  with  the 
maximum  temperature  T'f%.  The  depths  of  the  first  and 


d2F 

d2F 

d2F 

dz 2 

Z 1 

Z2 

dz 2 

-  dz2 

ZN 

Table  2.  Number  of  continental  shelf  salinity  profiles  in  Aug. 


Area 

DM  type 

SM  type 

Total 

SC  shelf 

205 

61 

266 

NC/WB  shelf 

87 

363 

450 

EB  shelf 

9 

100 

109 

Total 

311 

524 

825 

second  maxima  of  d2T/dz2  are  taken  as  top  d [  and  bot¬ 
tom  d 2  of  the  upper  and  lower  thermoclines  (Fig.  5c). 

The  SM-type  S  profiles  are  featured  by  a  maximum 
(minimum)  of  d2S/dz2  at  top  (bottom)  of  the  halocline. 
From  the  ocean  surface  downward,  the  depths  of  the 
first  maximum  and  the  first  minimum  of  d2S/dz2  are 
taken  as  top  df  and  bottom  d\  of  the  halocline  for  SM- 
type  S  profiles  (Fig.  5d). 


c.  Features  of  the  deep  water  profiles 

There  are  three  types  of  T  profiles  (DM,  SM,  and  AD 
types)  and  only  one  type  of  S  profile  (SM  type)  in  the 
deep  water  regime. 

The  DM-type  T  profiles  reveal  a  thick  surface  mixed 
layer  (Fig.  6a)  above  a  positive  gradient  thermocline. 
The  SM-type  '/’-profiles  are  characterized  by  a  minimum 
(maximum)  of  d2T/dz2  at  top  (bottom)  of  the  upper  ther¬ 
mocline.  Thus,  the  depths  of  the  first  minimum  and  the 
first  maximum  of  d2T/dz 2  are  taken  as  top  df  and  bottom 
d\  of  the  upper  thermocline  (Fig.  6b).  Below  the  upper 
thermocline,  a  region  of  cold,  nearly  isothermal  water 
exists,  which  is  a  remnant  of  the  previous  water’s  cool¬ 
ing  and  convective  mixing.  Beneath  this  isothermal  lay¬ 
er,  is  a  positive  gradient  lower  thermocline  similar  to 
that  of  the  DM  type. 

The  AD-type  T  profiles  are  depicted  by  maxima  of 
d2T/dz2  at  both  top  and  bottom  of  the  upper  and  lower 
thermoclines  and  a  depth,  df^,  with  a  maximum  tem¬ 
perature,  r^ax.  The  depths  of  the  first  and  second  max¬ 
ima  of  d2T/dz2  are  taken  as  top  df  and  bottom  df  of  the 
upper  thermocline  (Fig.  6c).  Below  the  depth  of  intru¬ 
sion  of  warm  Bering  Sea  water,  the  temperature  profile 
is  similar  to  that  described  by  the  DM  profiles. 

The  SM-type  S  profile  is  characterized  by  a  maximum 
(minimum)  of  d2S/dz2  at  top  (bottom)  of  the  upper  halo¬ 
cline.  From  the  ocean  surface  downward,  the  depths  of 
the  first  maximum  and  the  first  minimum  of  d2S/dz2  are 
taken  as  the  top  df  and  bottom  df  of  the  upper  halocline 
(Fig.  6d).  A  lower  halocline,  when  present,  is  often  con¬ 
current  with  the  positive  thermocline  of  the  SM-type  T 
profile  and  is  characterized  by  a  minimum  of  d2S/dz2  at 
the  lower  halocline. 


Table  1.  Number  of  continental  shelf  temperature  profiles  in  Aug. 


Area 

DM  type 

SM  type 

AD  type 

Total 

SC  shelf 

93 

228 

14 

335 

NC/WB  shelf 

140 

228 

174 

542 

EB  shelf 

2 

128 

1 

131 

Total 

235 

584 

189 

1008 

Table  3.  Seasonal  variation  of  number  of  deep  water  temperature 
profiles. 


Months 

DM  type 

SM  type 

AD  type 

Total 

Nov-Apr 

23 

0 

3 

26 

Aug 

28 

24 

10 

62 
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d.  Surface  mixed  layer 

The  bottom  of  the  mixed  layer  is  approximately  the 
same  as  the  top  of  the  upper  thermocline  (or  halocline). 
The  mixed  layer  has  a  depth  of  df  that  is, 

MLD'7  =  df. 

When  df  tends  to  an  infinitesimally  small  depth  (df 
—  0),  that  is,  the  upper  thermocline  (halocline)  reaches 
the  surface,  the  surface  mixed  layer  disappears.  We  call 
such  profiles  stratified  and  treat  them  as  a  special  case 
of  the  ordinary  SM-type  or  AD-type  profiles  when  df 
is  extremely  small. 

e.  Lower  thermocline 


of  the  entrainment  zone  equals  the  temperature  (sa¬ 
linity)  at  the  top  of  the  thermocline.  Therefore,  the 
temperature  and  salinity  jump  across  the  entrainment 
zone  is  represented  as 

A F  =  F( 0)  -  F(~df),  F  =  T,  S.  (7) 


c.  Upper  thermocline  and  halocline 

The  upper  thermocline  and  halocline  are  located  be¬ 
tween  z  =  ~ df  and  z  =  —df  Vertical  gradients  of  the 
upper  thermocline  (Gf)  and  halocline  (Gf)  for  SM-type 
T,  S  profiles  are 


The  lower  thermocline  appears  in  the  Beaufort  Sea 
deep  water.  The  major  geometric  feature  of  T  profiles 
is  the  existence  of  minimum  and  maximum  temperatures 
below  the  upper  thermocline  for  the  SM-type  and  AD- 
type  profiles  (Figs.  6b, c), 

Tmin  =  min  T{z)  and  7m.lx  =  max  T(z),  (3) 

z<-dj  z<~d\ 

dl  =  -\z\T{z)  =  TmJ,  (4) 

dl=  ~{z\T(z)  =  U,  (5) 

where  d\  and  d\  are  the  depths  of  the  minimum  and 
maximum  temperatures  below  the  upper  thermocline 
and  the  corresponding  data  are  T(—df)  and  T(—df).  We 
may  identify  the  layer  between  d\  and  d\  as  the  lower 
thermocline.  Usually,  between  the  upper  and  lower  ther- 
moclines  there  is  a  transition  zone  with  a  very  weak 
vertical  gradient. 

5.  Thermohaline  geometric  model 

Let  T  and  S  be  represented  by  F.  After  the  four  depths 
df  df  df  and  df  are  determined  from  temperature  and 
salinity  profiles,  we  may  obtain  some  important  physical 
parameters  from  each  profile. 


Gu  = 


Gu  = 

°ha 


T(-dl)  -  T(-dl) 
dT2  -  dl 

S(-df)  -  S(-dl ) 
dl  —  d  f 


and 


(8) 


The  vertical  gradient  of  the  upper  thermocline  (G'f)  for 
AD-type  T  profiles  is  estimated  by 


Gu  =  - 


i  n-di)  -  n-di  „)j + 


d] 


(9) 


d.  Lower  thermocline  and  halocline 

The  lower  thermocline  and  halocline  are  located  be¬ 
tween  z  =  ~dl  and  z  =  ~dl-  Vertical  gradients  of  the 
lower  thermocline  (Gt'h)  and  halocline  (GU  are  esti¬ 
mated  by 


Glh  = 


T(-dl)  -  T(-dl) 


and 


a.  Mixed  layer 

The  water  above  the  depth  of  d[  is  taken  as  the  mixed 
layer.  Within  the  mixed  layer,  the  temperature  (salinity) 
is  assumed  to  be  uniform  with  depth;  that  is,  the  mixed 
layer  temperature  is  the  same  as  the  sea  surface  tem¬ 
perature: 

F(z)  =  /TO),  -d{  <  z  <  0  (mixed  layer).  (6) 


b.  Entrainment  zone 

The  entrainment  zone  below  the  mixed  layer  is  as¬ 
sumed  to  be  infinitesimally  thin  near  the  depth  of 
df  The  temperature  (salinity)  at  the  top  of  the  en¬ 
trainment  zone  is  the  same  as  the  mixed  layer  tem¬ 
perature,  and  the  temperature  (salinity)  at  the  bottom 


G 


L 

ha 


S(-dj)  -  S(-dj) 
dl  —  dl 


(10) 


6.  Model  evaluation 

a.  Model  profile 

The  thermohaline  parametric  model  turns  any  profile 
into  physical  parameters:  four  depths  (Tf ,  dl,  dl,  df), 
temperature  or  salinity  jump  A F,  and  upper  and  lower 
thermocline  (halocline)  strength.  Taking  SM-type  T,  S 
profiles  as  an  example  (remember  that  the  DM  type  is 
a  special  case  of  the  SM  type),  we  construct  a  modeled 
profile  F(Zj),  to  evaluate  the  performance  of  the  ther¬ 
mohaline  parametric  model: 
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and 


Hz)  =  F(  0), 

F(z)  =  F(-d[)  -  *  +  d\Fm-di)  -  F(-df)], 

Hz)  =  Fi-do  -  fF+  rffF[*Wf)  -  *WQ]. 
%)  =  fwo  -  fF+  -  /'■(-rff)], 

J4  -  off 


~d(  <  z  <  0, 
-(/?  <  z  <  -df, 

-<if  <  z  <  -  cff, 

-(/£  <  z  <  -(if, 


F(z)  =  i7(-^) 


z  +  (if 
(if  -  (if 


[F(-(if)  -  F(— (if)], 


-(if  <  z  <  -(if, 


(11) 


which  is  vertically  uniform  in  the  mixed  layer  and  piece- 
wise  linear  with  depth  below  the  mixed  layer  to  the 
deepest  depth  of  the  observational  point  z  =  —(if.  No¬ 
tice  that  the  “model  profile”  takes  the  value  of 
F(—d()  at  the  mixed  layer  base  z  =  —(if. 

AD-type  profiles  can  also  be  represented  by  (11)  ex¬ 
cept  for  the  upper  thermocline  (—(if  Sz<  —df),  which 
should  be  parameterized  by 

Hz)  =  T(-dn  -  J  +  d[  \T(-d'U  -  T(-df)], 

«max  -  d{ 

-(/Lx  <  z  <  -(if, 


and 

Hz)  =  n-du  -  z  +  w(-di) 

dl  ~  “max 

-(if  <  7  <  —dT 
2  —  —  “max’ 


where  T(~dhmJ  = 


n-dL j], 
(12) 


b.  Root-mean-square  error 

The  model  performance  can  be  evaluated  by  the  root- 
mean-square  error  (rmse), 


rmse,„  =  ^  2  [F„,fc)  -  F„,(z;)]2,  (13) 

where  the  subscript  m  indicates  the  mth  profile  Fm(z). 
When  we  consider  a  group  of  profiles,  we  may  compute 
the  mean  value. 


rmse 


1 

M 


2 


rmse,„, 


(14) 


c.  Statistical  tests 

For  any  profile,  F(z,),  =  1,  2,  .  .  . ,  N,  using  (11), 
and/or  (11)  and  (12),  we  obtain  a  corresponding  model 
profile.  Fill),  i  =  1,  2,  .  .  .  ,  N.  Both  F(z)  and  Hz)  have 
the  same  values  at  the  depths  0,  —(if,  —(if,  —(if,  —(if, 
and  —(if,  but  do  not  necessarily  have  the  same  values 
at  the  other  observational  depths.  If  F(z)  fits  F(z.)  well, 
the  difference  SFiz.,)  =  Fiz.,)  —  /•(z,)  should  be  very 
small  at  all  observational  points.  We  need  to  test  if  the 
difference  is  small  enough  to  be  neglected.  The  first  test 
(paired  data  test)  is  used  to  justify  whether  the  mean 
value  of  SF(zi)  is  taken  as  zero.  The  second  test  (cor¬ 
relation  test)  is  utilized  to  see  whether  F(z)  correlates 
well  with  F(z) . 


1)  Paired  data  test  (r  test) 


After  calculating  the  mean  and  standard  deviation  of 
8F(z,), 


8F  =  hi  SF(zt), 

F V  i=  1 


S 


2 

SF 


1 

(N-  1) 


2 


[SF(Zl)  -  8FY. 


(15) 


we  begin  with  the  null  hypothesis  that  8F  is  zero.  The 
significance  level  (a  =  0.001)  is  the  probability  that  the 
given  value  of 


8F 

sSF/VN 


(16) 


is  exceeded  purely  by  chance.  This  value  satisfies  the  t 
distribution  with  N  —  1  degrees  of  freedom.  If  Id  > 
taN.  i ,  we  reject  the  null  hypothesis,  and  the  modeled 
profile  F(z)  does  not  fit  the  observation  F(z)  and  should 
be  rejected. 


where  M  is  total  number  of  profiles.  If  the  model  is 
perfect,  Fm(zt)  =  FJz,,).  The  rmsem  should  be  zero. 
Usually,  the  model  is  not  perfect.  The  smaller  the  rmse, 
the  better  the  performance  of  the  model. 


2)  Correlation  test  (F-test) 

After  a  modeled  profile  F(z)  passes  the  paired  data 
test,  we  start  with  the  null  hypothesis  that  F(z)  does  not 
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Table  4.  Number  of  fitted  profiles  and  rmse  before  and  after  the  t- 
and  F- tests. 


Shelf  water 

Deep  water 

Summer 

Summer 

Winter 

KK  rmse 

KK  rmse 

KK  rmse 

T 


Before 

1027 

0.25 

62 

0.31 

26 

0.14 

After 

850 

0.18 

61 

0.3 

21 

0.12 

Before 

788 

0.24 

19 

0.22 

24 

0.17 

After 

580 

0.16 

19 

0.22 

18 

0.14 

correlate  with  the  observation  F(z).  The  significance 
level  ( a  =  0.05)  is  the  probability  that  the  given  value  of 

1  [%)  -  h2 

f  =  -N - - -  (17) 

2  [%)  -  HzdYKN  -  2) 

i=0 

is  not  exceeded  purely  by  chance.  This  value  satisfies 
the  F  distribution  with  ( 1 ,  A'  —  2)  degrees  of  freedom. 
If/  >  F„(l,  N  —  2),  we  reject  the  null  hypothesis  and 
conclude  that  the  modeled  profile  F(z)  does  correlate 
with  the  observation  F(z). 

7.  Model-data  comparison 

a.  Overall  performance 

Any  profile  that  passes  both  the  f-test  and  the  f -test 
is  termed  a  fitted  profile.  The  more  fitted  profiles  that 
exist,  the  more  representative  they  are  of  the  geometric 
model.  Taking  all  the  August  profiles  as  an  example, 
almost  all  the  deep  water  profiles  are  fitted  profiles,  but 
only  about  80%  of  the  shelf  profiles  are  fitted  profiles 
(Table  4).  After  the  two  tests,  the  rmse  of  the  shelf 
profiles  has  drastic  reductions:  0.25°-0.18°C  and  0.24- 
0.16  psu  (Table  4).  We  checked  those  shelf  profiles  that 
did  not  fit  the  model  (around  20%)  and  found  that  the 
most  of  them  had  a  very  poor  vertical  resolution.  It  was 
hard  for  the  geometric  model  to  fit  those  profiles. 

b.  Shelf  water 

During  winter  (November-April)  most  T,  S  profiles 
exhibit  a  single-layer  structure  (DM-type  profiles).  It 
therefore  was  not  necessary  to  use  the  geometric  model 
to  parameterize  those  profiles. 

During  summer  (August)  T,  S  profiles  have  a  mul¬ 
tilayer  structure.  Figure  8  shows  the  model  performance 
of  the  SM-type  T,  S  profiles  for  the  three  shelf  regions: 
(a)  southern  Chukchi  (SC)  Sea,  (b)  northern  Chukchi/ 
western  Beaufort  (NC/WB)  Sea,  and  (c)  eastern  Beau¬ 
fort  (EB)  Sea.  After  application  of  the  t-  and  F-tests, 
the  mean  rmse  of  temperature  is  0.14°C  for  SC,  0.1 9°C 
for  NC/WB,  and  0.29°C  for  EB  (Table  5).  The  mean 


rmse  of  salinity  is  0.07  psu  for  SC,  0.20  psu  for  NC/ 
WB,  and  0.29  psu  for  EB  (Table  6).  Figure  9  shows  the 
model  performance  of  the  AD-type  T  profiles  for  the 
NC/WB  shelf  region.  After  applying  the  t-  and  F-tests, 
the  mean  rmse  of  temperature  is  0.3°C  for  NC/WB  (Ta¬ 
ble  5). 

c.  Deep  water 

The  deep  water  exhibits  a  multilayer  structure  all  year 
round.  Comparison  between  the  observed  and  modeled 
profiles  (Fig.  10)  shows  that  the  geometric  model  pa¬ 
rameterizes  the  various  type  profiles  reasonably  well. 
The  rmse  of  the  winter  T  profiles  is  0.1 2°C.  The  rmse 
of  the  summer  T  profiles  is  0.25°C  for  the  DM-type 
profiles,  0.28°C  for  the  SM-type  profiles,  and  0.45°C 
for  the  AD-type  profiles  (Table  5).  The  rmse  of  the  S 
profiles  is  0.14  psu  for  the  winter  profiles  and  0.22  for 
the  summer  profiles  (Table  6). 

The  upper  portion  of  the  Beaufort/Chukchi  Sea  deep 
water  experiences  a  seasonal  variation.  During  the  win¬ 
ter  season  (November-April),  the  surface  temperature 
mixed  layer  is  much  deeper  than  the  surface  salinity 
mixed  layer  (Fig.  10):  MLDtr)  =  170  m,  MLDIS>  =  30 
m.  During  other  seasons  (May-October),  the  surface 
mixed  layer  may  shoal.  August  T  profiles  show  the  ex¬ 
istence  of  two  types  of  thermohaline  structures:  (a)  deep 
mixing  (Fig.  lOe)  and  (b)  shallow  mixing  (Fig.  lOg)  and 
only  one  type  (shallow  mixing)  for  the  salinity  (Fig. 
lOi). 

8.  Beaufort/Chukchi  Sea  thermohaline  features 

The  thermohaline  geometric  model  turns  each  tem¬ 
perature  (or  salinity)  profile  into  a  set  of  physical  pa¬ 
rameters.  During  the  winter  the  water  throughout  the 
Beaufort/Chukchi  Sea  shelf  regions  is  quite  uniform  in 
the  vertical.  As  the  summer  approaches,  an  increase  in 
the  solar  insolation  begins  the  stratification  process  and 
causes  a  more  complicated  multilayer  thermohaline 
structure  in  the  Beaufort/Chukchi  Sea  shelf.  Therefore, 
we  choose  August  profiles  as  examples  for  illustrating 
the  usefulness  of  the  geometric  model. 

a.  SST  and  SSS 

Seasonal  evolution  of  SST  (Fig.  1 1)  constructed  from 
the  MOODS  data  (1970-93),  shows  some  interesting 
features.  Due  to  data  sparseness,  we  show  only  the  hor¬ 
izontal  distribution  of  warm  water  (SST  >  0°C)  and 
cold  water  (SST  <  0°C).  During  the  winter  season  (No¬ 
vember-April),  radiative  cooling  makes  the  whole 
Beaufort/Chukchi  SST  fall  below  0°C  (Fig.  1  la).  During 
the  spring  season  (May-July)  warm  water  appears  in 
the  southern  Chukchi  Sea  (Fig.  1  lb),  this  is  caused  by 
gradual  increase  of  the  solar  radiation  and  the  advected 
warm  water  through  the  eastern  Bering  Strait.  During 
summer  season  (August),  SST  of  the  deep  water  is  still 
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Fig.  8.  August  SM-type  observed  T,  S  profiles  and  corresponding  modeled  profiles  for  the  shelf  water:  (a)  observed  SC  temperature,  (b) 
modeled  SC  temperature,  (c)  observed  NC/WB  temperature,  (d)  modeled  NC/WB  temperature,  (e)  observed  EB  temperature,  (f)  modeled 
EB  temperature,  (g)  observed  SC  salinity,  (h)  modeled  SC  salinity,  (i)  observed  NC/WB  salinity,  (j)  modeled  NC/WB  salinity,  (k)  observed 
EB  salinity,  and  (1)  modeled  EB  salinity. 


cold  (below  0°C),  and  still  quite  cold  in  the  northern 
Chukchi/western  Beaufort  Sea  (Figs.  1  lc,d,e).  However, 
SST  of  the  shelf  water  is  warm  (above  0°C)  in  both 
southern  Chukchi  Sea  and  eastern  Beaufort  Sea.  During 
the  fall  season  (October),  the  cold  SST  almost  occupies 
the  whole  Beaufort  Sea.  However,  the  southern  Chukchi 
Sea  is  still  warm  (Fig.  Ilf). 

Similar  to  the  SST  fields,  we  show  the  seasonal  and 
horizontal  variations  of  salty  and  freshwater  (Fig.  12). 
The  criterion  for  separating  the  salty  and  freshwater  is 
31.5  psu  during  November-June,  and  26  psu  during 
July-October. 

Figures  11  and  12  clearly  show  the  intrusion  of  the 
warm  salty  Bering  Sea  water  into  the  Chukchi  Sea  and 


Table  5.  The  rmse  of  temperature  profiles  (°C). 


DM 

DM 

SM 

AD 

Area 

(winter) 

(summer) 

(summer) 

(summer) 

SC  shelf 

— 

— 

0.14 

— 

NC/WB  shelf 

— 

— 

0.19 

0.30 

EB  shelf 

— 

— 

0.29 

— 

Deep  water 

0.12 

0.25 

0.28 

0.45 

Beaufort  Sea  especially  from  May  to  October.  This  in¬ 
flux  was  discussed  by  Mountain  et  al.  (1976),  Paquette 
and  Bourke  (1974),  and  others.  Besides,  our  results  also 
show  strong  mixing  of  the  warm  salty  Bering  Sea  water 
mass  with  the  cold  fresh  Beaufort  Sea  water  during 
winter  season  (November-April). 

b.  Summer  (August)  shelf  water  thermohaline 
structure 

1)  Southern  Chukchi  Sea  Shelf 

Figure  13  exhibits  the  histograms  of  the  subsurface 


parameters  [MLD7 , 

<1 

-  d D,  Gth, 

MLD5,  AS, 

Table  6.  The 

rmse  of  salinity  profiles  (psu). 

DM 

DM 

SM 

Area 

(winter) 

(summer) 

(summer) 

SC  shelf 

— 

— 

0.07 

NC/WB  shelf 

— 

— 

0.20 

EB  shelf 

— 

— 

0.29 

Deep  water 

0.14 

— 

0.22 
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Fig.  9.  August  AD-type  observed  T  profiles  and  corresponding  modeled  profiles:  (a)  NC/WB  shelf  region  observed  profiles,  (b)  NC/WB 
shelf  region  modeled  profiles,  (c)  deep  water  observed  profiles,  and  (e)  deep  water  modeled  profiles.  The  rmse’s  are  listed  in  Table  5. 


(ds2  —  df),  Gha]  for  the  southern  Chukchi  Sea  shelf.  The 
horizontal  axis  in  each  figure  represents  the  correspond¬ 
ing  parameters.  The  vertical  axis  is  the  frequency.  There 
are  187  temperature  profiles  passing  both  t-  and  F-tests. 
The  thermal  mixed  layer  depth  (MLDr)  varies  from  0 
(no  thermal  mixed  layer)  to  40  m  with  a  mean  value  of 
10.1  m  (Fig.  13a),  the  temperature  jump  at  the  mixed 
layer  base  (AT)  varies  from  0°  (no  thermal  mixed  layer) 


to  4.5°C  with  a  mean  value  of  0.72°C  (Fig.  13b),  the 
thermocline  thickness  (d\  —  d[)  varies  from  2  to  44  m 
with  a  mean  value  of  11.9  m  (Fig.  13c),  and  the  ther¬ 
mocline  gradient  (Gth)  varies  from  near  0  to  1.75°C  m  1 
(Fig.  13d)  with  a  mean  value  of  0.43°  C  m_1. 

In  August,  there  are  205  DM-type  and  61  SM-type 
salinity  profiles  (Table  2).  For  the  SM-type  profiles, 
MLD5  changes  from  0  (no  salinity  mixed  layer)  to  14 
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Fig.  10.  Winter  (Nov- Apr)  and  summer  (Aug)  deep  water  observed  and  modeled  T,  S  profiles:  (a)  observed  winter  T  profiles,  (b)  modeled 
winter  T  profiles,  (c)  observed  winter  S  profiles,  (d)  modeled  winter  S  profiles,  (e)  observed  summer  DM-type  T  profiles,  (f )  modeled  summer 
DM-type  T  profiles,  (g)  observed  summer  SM-type  T  profiles,  (h)  modeled  summer  SM-type  T  profiles,  (i)  observed  summer  S  profiles,  and 
(j)  modeled  summer  S  profiles.  The  rmse’s  are  listed  in  Tables  5  and  6. 


m  with  a  mean  value  of  5.0  m  (Fig.  13e),  the  salinity 
jump  at  the  mixed  layer  base  (AS)  varies  from  0  (no 
salinity  mixed  layer)  to  1.8  psu  with  a  mean  value  of 
0.20  psu  (Fig.  13f),  the  halocline  thickness  ( d\  —  df) 
varies  from  3  to  45  m  with  a  mean  value  of  10.2  m 
(Fig.  13g),  and  the  halocline  gradient  (Glla)  varies  from 
near  0  to  1.48  psu  nu1  with  a  mean  value  of  0.40  psu 
m-1  (Fig.  13h). 

2)  Northern  Chukchi/western  Beaufort  Sea 
Shelf 

Figure  14  exhibits  the  thermohaline  features  for  the 
northern  Chukchi/western  Beaufort  Sea  shelf.  There  are 
162  temperature  profiles  passing  the  t-  and  E-tests. 
MLD7  varies  from  0  (no  thermal  mixed  layer)  to  28  m 
with  a  mean  value  of  7.7  m  (Fig.  14a),  AT  varies  from 
0°  (no  thermal  mixed  layer)  to  6.5°C  with  a  mean  value 
of  0.62°C  (Fig.  14b),  the  thermocline  thickness  (df  — 
df)  varies  from  3  to  47  m  with  a  mean  value  of  14.0  m 
(Fig.  14c),  and  Gth  varies  from  near  0  to  1.7°C  iru1  (Fig. 
14d)  with  a  mean  value  of  0.33°C  nu1.  There  are  252 


salinity  profiles  passing  the  t-  and  /--tests.  MLDS  chang¬ 
es  from  0  (no  salinity  mixed  layer)  to  20  m  with  a  mean 
value  of  4.7  m  (Fig.  14e),  AS  varies  from  0  (no  salinity 
mixed  layer)  to  7.8  psu  with  a  mean  value  of  0.55  psu 
(Fig.  14f),  the  halocline  thickness  (df  —  df)  varies  from 
2  to  110  m  with  a  mean  value  of  18.4  m  (Fig.  14g), 
and  Gha  varies  from  near  0  to  1.32  psu  iru1  with  a  mean 
value  of  0.36  psu  irr1  (Fig.  14h). 

3)  Eastern  Beaufort  Sea  Shelf 

Figure  1 5  exhibits  the  subsurface  parameters  for  the 
eastern  Beaufort  Sea  shelf.  There  are  128  temperature 
profiles  passing  the  t-  and  E-tests.  For  the  thermal  struc¬ 
ture,  MLDr  varies  from  0  (no  thermal  mixed  layer)  to 
15  m  with  a  mean  value  of  5.6  m  (Fig.  15a),  AT  varies 
from  0°  (no  thermal  mixed  layer)  to  4.4°C  with  a  mean 
value  of  0.43°C  (Fig.  15b),  the  thermocline  thickness 
(df  —  df)  varies  from  2  to  100  m  with  a  mean  value  of 
11.3  m  (Fig.  15c),  and  Gth  varies  from  near  0  to  1.83°C 
iru1  (Fig.  15d)  with  a  mean  value  of  0.65°C  nu1.  There 
are  100  salinity  profiles  passing  the  t-  and  E-tests.  For 
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Fig.  11.  Seasonal  variation  of  SST  presented  by  warm  water  (SST  >  0°C,  denoted  by  “•”)  and  cold  water  (SST  £  0°C,  denoted  by 

“+”).  Data  were  obtained  from  the  MOODS  for  1970-93. 


628 


JOURNAL  OF  ATMOSPHERIC  AND  OCEANIC  TECHNOLOGY 


Volume  16 


sBC1 1  -4(1 76/"+"sss<=31 .5o/oo) 


sBC5-6(26/"+"sss<=31 .5o/oo) 


-150 


-140 


Longitude(E),  Depth=-150(m) 


sBC-7(628/''o"sss<=26.o/oo) 


sBC-8(864/"o"sss<=26.o/oo) 


40 

Longitude(E),  Depth=-150(m) 


Longitude(E),  Depth=-150(m)  Longitude(E),  Depth=-150{m) 
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is  31.5  psu  for  Nov-Jun  and  26  psu  for  Jul-Oct.  Data  were  obtained  from  the  MOODS  for  the  years  1970-93. 


the  salinity  structure,  MLD  5  changes  from  0  (no  salinity 
mixed  layer)  to  12  m  with  a  mean  value  of  4.5  m  (Fig. 
15e),  AS  varies  from  0  (no  salinity  mixed  layer)  to  4.0 
psu  with  a  mean  value  of  0.40  psu  (Fig.  15f),  the  halo- 


cline  thickness  (d2  —  df)  changes  from  3  to  78  m  with 
a  mean  value  of  9.6  m  (Fig.  15g),  and  Gha  changes  from 
near  0  to  3.30  psu  m~'  with  a  mean  value  of  0.70  psu 
m-1  (Fig.  15h). 
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gradient  (°C  rn  1 ) ;  (e)  salinity  mixed  layer  depth  (m),  MLDS;  (f)  salinity  jump  at  mixed  layer  base  (psu),  AS;  (g)  upper-halocline  thickness 
(m);  and  (h)  halocline  gradient  [psu  m  ']. 


c.  Summer  (August)  deep  water  thermohaline 
structure 

Figure  16  exhibits  the  subsurface  parameters  for  the 
Beaufort  Sea  deep  water.  There  are  only  24  SM-type 
temperature  profiles  passing  the  t-  and  /-’-tests.  For  the 
thermal  structure,  MLD'  varies  from  0  (no  thermal 
mixed  layer)  to  5  m  with  a  mean  value  of  1.3  m  (Fig. 
16a),  AT  varies  from  0°  (no  thermal  mixed  layer)  to 
4.0°C  with  a  mean  value  of  0.26°C  (Fig.  16b),  the  upper 
thermocline  thickness  ( d\  —  df)  varies  from  2  to  240 
m  with  a  mean  value  of  67.3  m  (Fig.  16c),  and  G£  varies 
from  near  0  to  0.4°C  m  1  (Fig.  14d)  with  a  mean  value 
of  0.19°C  m  L 

There  are  19  salinity  profiles  passing  the  t-  and  F- 
test.  For  the  salinity  structure,  MLDS  changes  from  0 
(no  salinity  mixed  layer)  to  5  m  with  a  mean  value  of 
0.84  m  (Fig.  16e),  AS  varies  from  0  (no  salinity  mixed 
layer)  to  0.1  psu  with  a  mean  value  of  0.0053  psu  (Fig. 
16f),  the  upper  halocline  thickness  (d\  —  df)  varies  from 
9  to  25  m  with  a  mean  value  of  13.5  m  (Fig.  16g),  and 


varies  from  near  0  to  0.58  psu  m_1  with  a  mean 
value  of  0.31  psu  ur1  (Fig.  16h). 

These  thermohaline  features  of  both  shelf  water  and 
deep  water  are  outlined  in  Tables  7  and  8. 

9.  Conclusions 

Based  on  geometric  characteristics  of  profiles  (min¬ 
imum  and  maximum  d2T/dz2,  d2S/dz2),  several  key 
depths  such  as  MLD,  top  and  bottom  of  upper  ther¬ 
mocline,  and  halocline  can  be  identified.  Below  the  up¬ 
per  thermocline,  we  use  Tmin,  Tmax  to  determine  the  top 
and  bottom  of  the  lower  thermocline.  Other  thermo¬ 
haline  features  such  as  temperature  and  salinity  jumps 
at  the  mixed  layer  base  (AT,  AS),  upper-  and  lower- 
thermocline/halocline  gradients  can  also  be  identified. 

The  thermohaline  geometric  model  depicted  in  this 
paper  demonstrates  a  good  capability  to  compress  a 
large  dataset  into  a  small  enough  volume  of  physically 
meaningful  data  such  that  the  investigators  may  readily 
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Fig.  14.  Histograms  of  various  parameters  of  Aug  northern  Chukchi/western  Beaufort  Sea  shelf  thermohaline  feature  obtains  by  the 
geometric  model:  (a)  temperature  mixed  layer  depth  (m),  MI.D/  ;  (b)  temperature  jump  at  mixed  layer  base  (°C),  A T;  (c)  thermocline  thickness 
(m);  (d)  thermocline  gradient  (°C  m  1 ) ;  (e)  salinity  mixed  layer  depth  (m),  MLDS;  (f)  salinity  jump  at  mixed  layer  base  (psu),  AS;  (g) 
halocline  thickness  (m);  and  (h)  halocline  gradient  [psu  m  ']. 


assimilate  and  interpret  it.  Furthermore,  the  geometric 
model  can  process  multiple  type  profiles.  The  overall 
root-mean-square  errors  of  the  geometric  temperature 
model  is  0.18°C  for  the  shelf  water,  and  0.12°C  (winter) 
to  0.3°C  (summer)  for  the  deep  water.  The  overall  root- 
mean-square  error  of  the  geometric  salinity  model  is 
0.16  psu  for  the  shelf  water,  0.14  psu  (winter)  to  0.22 
psu  (summer)  for  the  deep  water. 

The  Beaufort/Chukchi  Sea  shelf  water  experiences  a 
strong  seasonal  variation.  The  geometric  model  pro¬ 
posed  here  clearly  shows  large  thermohaline  variabili¬ 
ties  and  the  multiple  structures.  During  the  winter  season 
(November-April),  the  surface  mixed  layer  for  both  T 
and  S  extends  from  the  surface  to  bottom.  During  the 
other  seasons  (May-October),  the  ocean  mixed  layer 
shoals.  Two  distinguishing  processes  (surface  warming 
and  horizontal  advection)  make  two  types  of  multilayer 
temperature  profiles  (shallow-mixing  type  caused  by  the 
surface  warming,  and  advection  type  caused  by  a  strong 
sublayer  advection)  and  only  one  type  of  multilayer  sa¬ 


linity  profiles  (shallow  mixing).  During  summer  (Au¬ 
gust),  the  mean  thermal  mixed  layer  depth  shoals  from 
10.1  m  for  the  southern  Chukchi  Sea  shelf  to  7.7  m  for 
the  northern  Chukchi/western  Beaufort  Sea  shelf,  and 
further  reduces  to  5.6  m  for  the  eastern  Beaufort  Sea 
shelf.  The  salinity  mixed  layer  depth  is  usually  thinner 
than  the  thermal  mixed  layer  depth  in  the  northern  Chuk¬ 
chi/western  Beaufort  and  the  eastern  Beaufort  Sea  shelf 
(mean  value  4.6  m),  but  thicker  than  the  thermal  mixed 
layer  in  the  southern  Chukchi  Sea,  where  77%  of  S 
profiles  have  a  one-layer  structure  (mixed  layer  depth 
is  the  same  as  the  water  depth). 

The  upper  portion  ( z  >  —150  m)  of  the  Beaufort/ 
Chukchi  Sea  deep  water  experiences  a  seasonal  varia¬ 
tion;  however,  the  lower  portion  (z  —  —150  m)  does 
not  have  an  evident  seasonal  variation.  During  the  win¬ 
ter  season  (November-April),  the  surface  thermal  mixed 
layer  (—170  m)  is  much  deeper  than  the  surface  salinity 
mixed  layer  (—30  m).  During  the  other  seasons  (May- 
October),  the  ocean  mixed  layer  shoals.  There  are  two 
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(d2mean=5.625  m) 


(tjumean=0.4252  C) 
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Fig.  15.  Histograms  of  various  parameters  of  Aug  eastern  Beaufort  Sea  shelf  thermohaline  feature  obtains  by  the  geometric  model:  (a) 
temperature  mixed  layer  depth  (m),  Mi.D7:  (b)  temperature  jump  at  mixed  layer  base  (°C),  A T:  (c)  thermocline  thickness  (m);  (d)  thermocline 
gradient  (°C  m  "■);  (e)  salinity  mixed  layer  depth  (m),  MLD5;  (f)  salinity  jump  at  mixed  layer  base  (psu),  AS;  (g)  halocline  thickness  (m); 
and  (h)  halocline  gradient  [psu  m  ']. 


types  of  multilayer  temperature  profiles  (shallow-mix¬ 
ing  and  advection  types)  and  only  one  type  of  multilayer 
salinity  profiles  (shallow  mixing).  During  summer  (Au¬ 
gust),  the  upper  thermocline  and  halocline  almost  reach 
the  surface  for  the  shallow  mixing  type  T,  S  profiles. 
The  mean  thermal  and  salinity  mixed  layer  depths  are 
1.3  and  0.8  m,  respectively. 

Strong  temperature  and  salinity  jumps  at  the  mixed 
layer  base  appear  in  the  summer  season.  The  mean  val¬ 
ues  of  AT,  AS  of  August  profiles  are  around  0.72°C, 
0.20  psu  in  the  southern  Chukchi  Sea  shelf;  0.62°C,  0.55 


psu  in  the  northern  Chukchi-western  Beaufort  Sea  shelf; 
and  0.43°C,  0.40  psu  in  the  eastern  Beaufort  Sea  shelf. 

Strong  upper  thermocline/halocline  gradients  are 
found  in  the  summer  season.  Both  upper-thermocline 
and  halocline  gradients  (Gj(,  G,'j)  weaken  from  the 
southern  Chukchi  Sea  shelf  [0.43°C  m  1 ,  0.40  psu  m~'] 
to  the  northern  Chukchi/western  Beaufort  Sea  shelf 
[0.33°C  m_1,  0.36  psu  m_1],  and  then  strengthen  from 
the  northern  Chukchi/western  Beaufort  Sea  shelf  to  the 
eastern  Beaufort  Sea  shelf  [0.65°C/m~\  0.70  psu  m~*]. 


Table  7.  Mean  values  of  thermal  parameters  for  three  shelf 
regions  and  deep  water  in  Aug. 


MLDr  (m) 

A T  (°C) 

di  -  d\ 

Gu 

v-7lh 

SC  shelf 

10.1 

0.72 

11.9 

0.43 

NC/WB  shelf 

7.7 

0.62 

14.0 

0.33 

EB  shelf 

5.6 

0.43 

11.3 

0.65 

Deep  water 

1.3 

0.26 

67.3 

0.19 

Table  8.  Mean  values  of  haline  parameters  for  three  shelf  regions 
and  deep  water  in  Aug. 


MLD"  (m) 

AS  (psu) 

ds2  -  d\ 

GS, 

(psu  m  *) 

SC  shelf 

5.1 

0.20 

10.2 

0.40 

NC/WB  shelf 

4.7 

0.55 

18.4 

0.36 

EB  shelf 

4.5 

0.40 

9.6 

0.70 

Deep  water 

0.8 

0.01 

13.5 

0.31 
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(d2mean=1 .25  m)  (tjumean=0.2583  c) 
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Fig.  16.  Histograms  of  various  parameters  of  Aug  Beaufort/Chukchi  Sea  deep  water  SM-type  profiles  obtains  by  the  geometric  model: 
(a)  temperature  mixed  layer  depth  (m),  MLD7";  (b)  temperature  jump  at  mixed  layer  base  (°C),  AT;  (c)  upper-thermocline  thickness  (m);  (d) 
upper-thermocline  gradient  (°C  m_l);  (e)  salinity  mixed  layer  depth  (m),  MLDS;  (f)  salinity  jump  at  mixed  layer  base  (psu),  AS;  (g)  upper- 
halocline  thickness  (m);  and  (h)  upper-halocline  gradient  [psu  m  ‘j. 


It  is  a  challenge  to  interpolate/extrapolate  profile  data 
horizontally  and  the  profile  parameterization  may  help. 
Thus,  the  geometric  (or  parametric)  model  is  useful  in 
building  temperature  and  salinity  climatology  such  as 
GDEM. 

This  geometric  model  can  also  be  used  to  process 
output  from  any  dynamical  model  and  to  obtain  mean¬ 
ingful  products,  that  is,  mixed  layer  depth,  upper-  and 
lower-thermocline/halocline  depths  and  strengths,  and 
deep  layer  stratification. 
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